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The matrix element of B° — B° mixing is evaluated from QCD sum rules for three-point correlators with 
next-to-leading order accuracy in pQCD 1 . 
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1. Introduction 

Particle-antiparticle mixing in systems of neu- 
tral mesons of different flavors 

sd: K°-K°; cu : D°-D°; bd : B° - B° 

is the primary source of CP violation studies. 
Historically, the study of K° — K° mixing was 
very fruitful for elementary particle physics. It 
provided deep insights into delicate questions of 
weak interactions and gave first convincing proof 
of possibility for CP violation. The thorough 
quantitative analysis of mixing in the kaon system 
strongly constrained the physics of heavy parti- 
cles and possible scenarios for the extension of 
the light quark sector of the theory. The nu- 
merical value of the mass splitting between the 
eigenstates of the effective Hamiltonian for the 
kaon pair has been used to estimate the numerical 
value of the charm quark mass from the require- 
ment of GIM cancellation before the experimental 
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discovery of charm (see e.g. 2 ). The proof of ex- 
isting CP violation in the interaction of particles 
happens to be very important for our understand- 
ing of the structure of Universe and its evolution. 
For example, the violation of CP symmetry is one 
of necessary conditions for generating the baryon 
asymmetry of the Universe [3] . 

Presently the experimental studies of phenom- 
ena of CP violation and mixing are active in the 
area of heavy mesons for which they are con- 
sidered more promising. The charmed mesons 
D(uc) have recently attracted some interest and 
provided encouraging results |3j. The systems 
of Bd(db) and B s (sb) mesons are the current 
laboratory for performing a precision analysis 
of CP violation and mixing both experimen- 
tally and theoretically joltij . The experimental 
study is now under way in dedicated experiments 
at SLAC (BABAR collaboration, e.g. [7fB] ) and 
KEK (BELLE collaboration e.g. pTflj V 

2. Phenomenology of B° — B° mixing 

Mixing in a system of neutral pseudoscalar 
mesons is phenomenologically described by a 
2x2 effective Hamiltonian H e ff = (M — 
{hj} = {1)2} where M is related to the mass 
spectrum of the system and T describes the 
widths of the mesons. The time evolution of the 
system state vector (B°,B°) is governed by the 
equation 

The AB = 2 flavor violating interactions gener- 
ate the non-diagonal terms in the effective Hamil- 
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tonian. The mass difference Am 



M, 



heavy 



ight ~ 2|M 12 | has been preciseiy measured 
Am = 0.489±0.005(stat)±0.007(sysi) ps" 1 [TT] . 
This is an important observable which can be 
used to extract the top quark CKM parameters 
provided that the theoretical formulae are suffi- 
ciently accurate. The main difficulty of the anal- 
ysis is to account for the effects of strong interac- 
tions. 

The effective AB = 2 Hamiltonian is known 
at next-to-leading order (NLO) in QCD pertur- 
bation theory of the Standard Model JTSJ 



H 



AB=2 

eff 
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00 
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w (v tb *v td y VB So{x t ) (i) 
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4tt 
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Here Gf is a Fermi constant, M\y is the VF-boson 
mass, rjB = 0.55 ± 0.1 [Ej, J 5 = 1.627 in the 
naive dimensional regularization (NDR) scheme, 
So(x t ) is the Inami-Lim function ^3], and O(fi) = 
ipLladL){bLladL){^) is a local four-quark oper- 
ator at the normalization point p. Mass splitting 
of heavy and light mass eigenstates is 



Am = 2|(B | J H^j? =2 |B°>|. 



(2) 



The largest uncertainty in the calculation of the 
mass splitting is introduced by the hadronic ma- 
trix element A = (B \O(ii)\B°) that is poorly 
known 

3. Hadronic Bb parameter 

Theoretical evaluation of the hadronic matrix 
element A is a genuine non-perturbative task 
which should be approached with some non-direct 
techniques. The simplest approach ("factoriza- 
tion" J^j) reduces the matrix element A to the 
product of simpler matrix elements measured in 
leptonic B decays 

A f = l(B o \b Lla d L \0)(0\b Ll °d L \B o ) = 2 -f B ml 

where the leptonic decay constant fg is defined 
by the relation (0|6 L 7 M d L |.B o (p)) = ip^/s/2 and 
tub is the B° meson mass. A deviation from the 
factorization ansatz is usually described by the 



parameter B B defined as A = BbA^; in factor- 
ization Bb = 1. The evaluation of this parame- 
ter (and the analogous parameter B^ of K° — K° 
mixing) has long history. Many different results 
were obtained within approaches based on quark 
models, unitarity, ChPT. The approach of direct 
numerical evaluation on the lattice has also been 
used. The corresponding results can be found in 
the literature |16ll7ll8ll()W)l2ll2^l2\^4| . 

In my talk I report on the results of the cal- 
culation of the hadronic mixing matrix elements 
using Operator Product Expansion (OPE) and 
QCD sum rule techniques for three-point func- 
tions 17 18 2 5l2f)l27j . This approach is very close 
in spirit to lattice computations |23| . The advan- 
tages of this approach are: 

• It is a model- independent, first-principles 
method. The difference with the lattice ap- 
proach is that the QCD sum rule approach 
uses an asymptotic expansions of a Green's 
function computed analytically while on the 
lattice the function itself can be numeri- 
cally computed provided the accuracy of the 
technique is sufficient. 

• The sum rule techniques provide a consis- 
tent way of taking into account perturba- 
tive corrections which is needed to restore 
the RG invariance of physical observables 
usually violated in the factorization approx- 
imation |25| . 

A concrete realization of the sum rule method 
used in the analysis consists in the calculation of 
the moments of the three-point correlation func- 
tion of the interpolating operators of the i?-meson 
and the local operator O(n) responsible for the 
B° — B° transitions. 

4. pQCD analysis at three loops 

Consider the three-point correlation function 

n( P i,p 2 ) (3) 

= J dxdy(0\TJ B (x)O(0)J B (y)\0)e^ x -^ 

of the relevant AB — 2 operator <D(fi) and inter- 
polating currents for the B -meson Jb = (m/, + 



Sum rules for B — B mixing 





Qg, 



Figure 1. Perturbation theory diagram at LO 



Figure 2. Factorizable diagrams at NLO 



mdjdi'ysb. Here m& is the b quark mass. The cur- 
rent J B is RG invariant and J B = ^(^7^756). 
The main relevant property of this current is 
(0|Js(0)|B°(p)> = f B m% where m B is the B- 
meson mass. A dispersive representation of the 
correlator reads 



n(pi,p 2 ) 



p{si,s 2l q 2 )dsids 2 

(*1-P?)(*2 



(4) 



where q = p 2 — p\. For the analysis of S — 
B° mixing within the sum rule framework this 
correlator can be computed at q 2 = 0. 

Phenomenologically the mixing matrix element 
determines the contribution of the i3-mesons in 
the form of a double pole to the three-point cor- 
relator 



n(pi,p 2 ) 



™ B -Pi 



Pi 



Because of technical difficulties of calculation, a 
practical way of extracting the B° — B° matrix 
element is to analyze the moments of the correla- 



tion function at pf 



p\ 



at the point q~ = 



M(i,j) = 



p(si,S2,0)dsids 2 

*1 b 2 



A theoretical computation of these moments re- 
duces to an evaluation of single-scale vacuum dia- 
grams and can be done analytically with available 
tools for the automatic computation of multi-loop 
diagrams. Note that masses of light quarks are 
small (e.g. 28 29I30| 1 and can be accounted for 
as small perturbation. This makes it possible to 
analyze the problem of mixing for B s mesons . 
The leading contribution to the asymptotic ex- 
pansion is given by the diagram shown in Fig. ^ 



At the leading order in QCD perturbation the- 
ory the three-point function T\(pi,p 2 ) of Eq. J3J 
completely factorizes 



n(pi,p2) = ^n M ( Pl )iP(p 2 ) 



(5) 



into a product of the two-point correlators IT^ (p) 

n»= P/ rV) (6) 

= / dxe ipx (0\TJ s (x)b L ^d L (0)\0). 



At the LO the calculation of moments is straight- 
forward since the double spectral density is ex- 
plicitly known in this approximation. Indeed, us- 
ing dispersion relation for the two-point correla- 



tor 



P(s)- 



30 p(s)ds 

,2 S — p 2 ' 



(7) 



16?r 2 



1 



nr 

s 



one obtains the LO double spectral density in a 
factorized form 

p LO (s 1 ,s 2 ,q 2 ) = -{pi ■P 2 )p(si)p(s 2 ). 

Thus, all PT contributions are of the factor- 
izable form at the leading order. First non- 
factorizable contributions to Eq. (@J appear at 
NLO. Of course, at NLO there are also the fac- 
torizable diagrams. Note that the classification of 
diagrams in terms of their factorizability is con- 
sistent as both classes are independently gauge 
and RG invariant. Consider first the NLO factor- 
izable contributions that are given by the product 
of two-point correlation functions from Eq. (jSJ, as 
shown in Fig. [21 Analytical expression for such 
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Figure 3. A non-factorizable diagram at NLO 



contributions can be obtained as follows. Writing 
n(p 2 ) = n LO (p 2 ) + n NLO (;p 2 ) one finds 



n NLo(Pl'P2) = g(Pl-P2)(n L o(Pl)IlNLo(P2) 

+n N Lo(Pi)n L o(P2))- 

The spectral density of the correlator IInlo(p 2 ) 
is known analytically that solves the problem 
of the NLO analysis in factorization. Even a 
NNLO analysis of factorizable diagrams is pos- 
sible as several moments of two-point correlators 
are known analytically |32| . 

The NLO analysis of non-factorizable contri- 
butions within perturbation theory is the main 
point of my talk. The analysis amounts to the 
calculation of a set of three-loop diagrams (a typ- 
ical diagram is presented in Fig. |3J. These di- 
agrams have been computed using the package 
MATAD for automatic calculation of Feynman 
diagrams [3*3"]. The package is applicable only for 
computation of scalar integrals, the decomposi- 
tion of the three-point amplitude into scalars is 
known [!""f]. The expression for the "theoretical" 
moments is 



M th (i,j) 



and represent 



where the quantities a^, b{j 
LO, NLO factorizable and NLO nonfactorizable 
contributions as shown in Figs. 11131 The NLO 
nonfactorizable contributions 6™ with i + j < 7 
are analytically calculated in ref. [I] for the first 
time. The calculation has been done with the sys- 
tem FORM [""""] and required about 24 hours of 
computing time on a dual-CPU 2 GHz Intel Xeon 



Figure 4. A model of phcnomenological spectrum 



machine. The analytical result for the lowest fi- 
nite moment M t ^(2,2) reads 



1 



022 



(16vr 2 ) 5 



b f - — 



40 16tt 2 



(9) 



f 8366187 
22 17500 



-C; 



84608 
875 



3 9 ' 
,33197 426319 



52500 315000 



Here S 2 = ^Cl 2 (|) = 0.2604..., ( 3 = 
C(3), and /i 2 = to 2 . Higher moments con- 
tain the same transcendental entries S2, C3i 71-2 
with different numerical coefficients. The numer- 
ical values for the moments are b™J : ^2(2345) = 
{0.68,1.22,1.44,1.56} and b^L, = {1.96,2.25}. 
The above theoretical results are used to extract 
the non-perturbative parameter Bb from the sum 
rules analysis. 



5. Results for Bb 

The "phcnomenological" side of the sum rules 
reads 



M ph {i,j) = -B B 



B 



(10) 



where the contribution of the B-meson is dis- 
played explicitly. The remaining parts are the 
contributions due to higher resonances and the 
continuum which are suppressed due to the mass 
gap A in the spectrum model. A rough picture of 
the phcnomenological spectrum is given in Fig. [3] 
For comparison we consider the factorizable ap- 
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proximation for both "theoretical" 



ML (i,j)= 1 + -^k 

th m 2(»+j) V 47T 



(11) 



and "phenomenological" moments, which, by 
construction, are built from the moments of the 
two-point function of Eq. © 



<(M)=o 



8 !%m B 



(12) 



The standard comparison of theoretical calcula- 
tion with the phenomenological representation al- 
lows to extract the residue in the senior pole and, 
therefore B B ■ Eq. (fTTfl and Eq. JHJ differ only due 
to non-factorizable corrections. 

To find the nonfactorizable addition to Bb 
from the sum rules we form ratios of the total 
and factorizable contributions. On the "theoreti- 
cal" side one finds 



M th (i,j) 



unf 



M t {(i,j) 



(13) 



This ratio is mass- independent. On the "phe- 
nomenological" side we have 



M ph (i,j) = B B + R B {zi +z i ) + C B z i +j 
M f ph {i,j) ~ 1 + R f {z J +z l ) + Cfz'+'J 



(14) 



where z — m B /(m B + A) is a parameter that 
describes the suppression of higher state contri- 
butions. A is a gap between the squared masses 
of the B-meson and higher states. R B , C B , R? 
and are parameters of the model for higher 
state contributions within the sum rule approach. 
In order to extract the non-factorizable contribu- 
tion to B B we write B B = 1 + AB. Similarly, one 
can parameterize contributions to "phenomeno- 
logical" moments due to higher f?-meson states 
by writing R B = R* + AR and C B = C f + AC. 
Clearly, AB = AR = AC = in factorization. 
The final formula for the determination of AB 
reads 



—6™/ = AB + AR(z J ~ 2 1 ~ 1 - 2 ^ 



+ z l 



ACz l+j - 4 



where AR and AC are free parameters of the fit. 
We take A = 0Am B for the B meson two-point 



correlator. This corresponds to the duality inter- 
val of 1 GeV in energy scale for the analysis based 
on finite energy sum rules |36| . The actual value 
of AB has been extracted using the least-square 
fit of all available moments. Estimating all uncer- 
tainties we finally find the NLO non-factorizable 
QCD corrections to AB due to perturbative con- 
tributions to the sum rules to be 



AB 



±1) 



a s (m) 

47T 



We checked the stability of the sum rules which 
lead to a prediction of AB. For mj = 4.8 GeV 
and a s (m b ) = 0.2 37 one finds AB = 0.1. 

6. Conclusion 

To conclude, the B° — B° mixing matrix el- 
ement has been evaluated in the framework of 
QCD sum rules for three-point functions at NLO 
in perturbative QCD. The effect of radiative cor- 
rections on is under complete control within 
pQCD and amounts to approximately +10% of 
the factorized value. The calculation can be fur- 
ther improved with the evaluation of higher mo- 
ments. The result is sensitive to the parameter z 
or to the magnitude of the mass gap A used in 
the parametrization of the spectrum. 
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